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ABSTRACT 

Muon catalyzed fusion is a process in which a negatively charged muon combines with two 
nuclei of isotopes of hydrogen, e.g. a proton and a deuteron or a deuteron and a triton, to form 
a muonic molecular ion in which the binding is so tight that nuclear fusion occurs. The muon 
is normally released after fusion has taken place and so can catalyze further fusions. As the 
muon has a mean lifetime of 2.2 microseconds, this is the maximum period over which a muon 
can participate in this process. This article gives an outline of the history of muon catalyzed 
fusion from 1947, when it was first realised that such a process might occur, to the present day. 
It includes a description of the contribution that Drachman has made to the theory of muon 
catalyzed fusion and the influence this has had on the author’s research. 

INTRODUCTION 

I am very grateful to the organizers for the invitation to speak at this Symposium to mark the 
retirement of Aaron Temkin and Dick Drachman. It is a very great pleasure to be able to 
contribute to this occasion in this way. 

The particle the muon that is central to muon catalyzed fusion was discovered by Anderson 
and his first graduate student, Neddermeyer, in 1936 when studying cosmic rays (ref. 1). In 
accordance with relativistic quantum field theory, there is a negatively charged muon, , and 
its positively charged antiparticle, /i + , each having mass 207m c , where m e is the mass of the 
electron. Each has a mean lifetime of 2.2 microseconds, before decaying through the weak 
interaction into an electron or a positron and a neutrino and an antineutrino. See, for example, 
ref. 2. The muon is a lepton, as it is unaffected by the strong interaction. It follows that it is a 
fermion with spin In what follows, I will assume that the muon being considered is 
negatively charged. 

If the electron in a hydrogen atom (H) is replaced by a muon, the result is a hydrogen-like 

atom p/x called muonium. The reduced mass of this atom is 186m e . Thus its bohr radius is 

7 “;, where a 0 is the bohr radius of H. The binding energy of each of its bound states is 186 x 
186 

(the corresponding value for the H atom). Muonium is thus very compact and strongly bound. 
The proximity of the positive charge on the proton and the negative charge on the muon make 
it similar to a neutron. 

An electron can bind two protons to form a weakly bound ion, Hj, the hydrogen molecular 
ion. If the electron in this ion is replaced by a muon, the resulting ion in its ground state is 
very compact and strongly bound like muonium. As is to be expected, these properties remain 
if either or both protons are replaced by a deuteron ( d ) or a triton (i). 

It has been known since the 1930s that the sun generates its energy by a fusion reaction in 
which hydrogen is converted into helium (see, for example, Bethe (ref. 3). However, this 
reaction requires a temperature of several million degrees to overcome the Coulombic repulsion 
between the protons involved in the reaction. The proton-proton fusion process is very much 
slower than those involving deuterons or tritons as it involves the weak interaction. Extensive 
research into using a fusion process of this type as a source of energy here on Earth has been 
going on for many years. However, even for the most favourable fusion process, d + i, very high 
temperatures are required. See, for example, ref. 4. This has made progress slow. 
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A natural question to ask is whether any way could be found of bringing about fusion that did 
not require very high temperatures. In 1947, Frank (ref. 5,6) suggested that in the presence of 
protons and deuterons a slow muon might bind to a p to form pp, which as noted earlier is 
similar to a neutron. This could come close to a d and form pdp which we have seen is tightly 
bound. Frank thought that it was sufficiently tightly bound that p-d fusion might occur. With 
luck, this would leave the muon free to catalyze further fusions, i.e. bring about p-d fusion 
while remaining unchanged at the end of the reaction. 

Shortly afterwards Sakharov (ref. 7) discussed the possibility of energy production by such a 
process and later it was considered further by Zel’dovich (ref. 8). Both concluded 
pessimistically that a muon was only likely to be able to catalyze one or two fusions. 

The first experimental observation of muon catalyzed fusion was made by Alvarez et al. in 1956 
(ref. 9). Though unaware of the above theoretical speculations, they correctly interpreted tracks 
in their bubble chamber at Berkeley as representing p-d fusion catalyzed by a single muon. 
Bubble chamber tracks representing the path of a muon in a mixture containing hydrogen and 
a small amount of deuterium that catalyzes two p-d fusions are shown in Figure 1 in ref. 10. 

Upon reading about Alvarez’s exciting discovery in the New York Times, Jackson (ref. 11) 
proceeded to make an analysis of energy production. Like Sakharov and Zel’dovich, he came to 
a pessimistic conclusion: very few fusions could be catalyzed by a single muon. 

The scene now shifted, after a pause, to the Soviet Joint Institute for Nuclear Research (JINR) 
at Dubna, north of Moscow. Interest in muon catalyzed fusion revived when experiments on 
d-d fusion catalyzed by a muon carried out by Dzhelepov et al. (ref, 12) at Dubna in 1966 
revealed a strong and unexpected temperature dependence in the rate of formation of ddp. 

This was a very exciting discovery. It strongly suggested that a resonant process was involved 
and held out the possibility of a large increase in the fusion rate under suitable conditions. 

It was not long before a form was suggested for this resonant process. In 1967, Vesman (ref. 13) 
proposed that the formation of ddp could occur by the reaction 


dp + D 2 — * [(ddp)*dee\ 


( 1 ) 


where the species on the right-hand side is a muonic molecular complex in which ddp , in a 
weakly bound excited state, forms one of the nuclei. The energy lost by the system through the 
formation of the weakly bound state goes into exciting the vibrational and rotational states of 
the muonic molecular complex. Due to the quantised structure of these states, this will be a 
resonant process and the formation rate will be very sensitive to the kinetic energy of the dp. 
This accounted qualitatively for the temperature dependence observed by Dzhelepov et al. 


Such a mechanism depends crucially on the existence of a weakly bound state with binding 
energy less than the 4.5 eV dissociation energy of D 2 . The problem was taken up by two 
Russian physicists, Gershtein and Ponomarev. Gershtein had earlier been involved in work on 
muon catalyzed fusion with ZePdovich (ref. 14). Ponomarev was his graduate student. 

After ten years’ work, Gershtein and Ponomarev and their group (refs. 15,16) were indeed able 
to confirm the existence of a weakly bound state of ddp . They showed that it had rotational 
and ‘vibrational’ quantum numbers ( J, v) = (1,1) and binding energy ~2 eV. 


The Born-Oppenheimer approximation does not work well for ddp as the ratio of the masses of 
the muon and the deuteron is too large. The method Ponomarev and his group used to remove 
this difficulty was the adiabatic representation method (ref. 17) in which the energy and wave 


function of the state were calculated using a large basis set of the form 




where 


r is the position vector of the muon with respect to the geometric center of the nuclei as origin, 
R is the internuclear vector, ^i(r) is a Born-Oppenheimer wave function for the muon and 
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Xi(R)/R is a nuclear wave function. The (xi(R)} were determined by solving close-coupling 
type integro-differential equations. 

Gershtein and Ponomarev were also able to show by this method that a corresponding weakly 
bound state of dtp exists with binding energy ~ 1 eV. In addition , they were able to show that 
rapid formation of dtfi by the Vesman mechanism would enable a muon to catalyze ^ 100 d-t 
fusions (ref. 15). The rapid formation rate of dtp was confirmed experimentally by Bystritsky 
et al. (ref. 18) in 1980. 

GROWTH OF INTEREST IN MUON CATALYZED FUSION 

Initially at Dubna and later at the Kurchatov Institute in Moscow, Ponomarev built up an 
extensive group of very able Russian physicists working on the problem of muon catalyzed 
fusion. The very promising advances described above took place at a time when the Cold War 
made direct contacts with Russian scientists difficult. However, the results were available in the 
literature, as can be seen from the list of references. They appeared in English after a slight 
delay, but greater intimacy could be achieved by a reader with a knowledge of Russian. 

In 1979 the accelerator at Dubna that had been used in the experiments was shut down for 
refurbishment. The experimental lead was taken up by collaborations at Los Alamos, led by 
Jones, at the Paul Scherrer Institute (PSI) in Switzerland led by Petitjean and at the 
University of Tokyo and KEK in Japan led by Nagamine. 

On the theoretical side, the work of Gershtein and Ponomarev became widely known. Outside 
the Soviet Union, expertise had accumulated for many years on calculations on few-body 
Coulombic systems using the Rayleigh-Ritz variational method. It was possible that this 
method might give more accurate results for the binding energies of states of ddp and dtp than 
the adiabatic representation method (ref. 17). 

The first such calculation in the West was carried out by Bhatia and Drachman (ref. 19). I 
consider that it was very much in character for Drachman to be the first in this way. I can 
think of a number of reasons why muon catalyzed fusion appealed to him. The Vesman 
mechanism (ref. 13), as substantiated by Gershtein and Ponomarev, was a very interesting and 
exciting process involving an unusual particle, a muon, and the formation of a muonic 
molecular complex associated with resonant behaviour. Muon catalyzed fusion was a novel 
topic to an extent that few topics are. It transcended more than national boundaries in the 
best spirit of science. Drachman’s knowledge of Russian made it possible for him to get 
up-to-date earlier on new developments in Soviet journals and go further and read items such 
as untranslated reports written by Ponomarev’s group. Finally, muon catalyzed fusion just 
might be a viable way of bringing about fusion at room temperature that could be used as a 
commercial source of energy. 

Bhatia and Drachman (ref. 19) carried out calculations of binding energies of all states of the 
muonic molecular ions ppp, pdp, ptp, ddp, dtp and ttp 7 except the states of ddp, dtp and ttp 
that have J > 1. They used the Rayleigh-Ritz variational method with up to 440 basis 
functions {^} containing Hylleraas-type functions, i.e. functions that are linear in one or more 
interparticle distances. 

For J — 0, 

A - e-^'+^r'rZVk D°+ (l, m, n > 0), 

where 

r a , 7*2 = distances of the muon from the nucleons, 

T \2 = distance between the nucleons, 
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and Dq + is a rotational harmonic (ref. 20). 

For J — 1, the rotational harmonics D] + and Dj“ and functions of the angle between the 
directions of v\ and r 2 were used to include the unit of angular momentum. 

In some cases, the calculated binding energy proved to be rapidly convergent and lower 
energies were obtained than in any previous calculations. For the key weakly bound states, 
however, they obtained binding energies consistent with previous results, but significantly 
smaller in magnitude. 

Bhatia and Drachman pointed out that similar results were obtained by Frolov and Efros 
(ref. 21) at the Kurchatov Institute in Moscow using up to 375 basis functions {xt} of the form 

Xt = exp[-Q U ri - a 2i r 2 - ai 2 ,fi 2 ]fi/{ri, r 2 ), 

where ft/ (r 1 , 1 * 2 ) is a function of r 1 and r 2 appropriate to the angular momentum J of the 
state under consideration. (This work was later published in ref. 22.) 

The calculations were followed by many subsequent calculations by the Rayleigh-Ritz 
variational method using various types of basis set. 

For example, Hu (ref. 23) and Szalewicz et al. (ref. 24) used Hylleraas-type basis functions 
similar to Bhatia and Drachman (ref. 19). Aissing and Monkhorst (ref. 25) relied on exponential 
basis functions similar to Frolov and Efros (refs 21,22) with many non-linear parameters. 
Kamimura (ref. 26) showed the effectiveness of using basis functions similar to Frolov and Efros 
but with exponentials in the form of Gaussians. Vinitsky et al. (ref. 27) and Hara and Ishihara 
(ref. 28) used Hylleraas-type basis functions with the muon part of the functions expressed in 
terms of prolate spheroidal coordinates. These calculations gave highly accurate results. 


THE RESONANT REACTION THAT LEADS TO FUSION 

The most favourable muonic molecular ion for muon catalyzed fusion is dtp. For example* the 
d-t fusion rate for this ion is more than 1000 times larger than for any other ion. Also, it 
produces nearly the highest energy per fusion (17.6 MeV). See, for example, ref. 29. 

In the resonant reaction 

tp + D 2 — > [(dtp)ndee] (2) 

dtp is formed in its J = 1, v = 1 state. However, rapid fusion only takes place if dtp is in a 
state with J — 0. It is thus important to know the binding energies and wave functions of the 
various states of dtp below the very weakly bound (1,1) state. These binding energies are 
shown in Figure 1 which is taken from a revue article on muon catalyzed fusion by Bhatia and 
Drachman (ref. 30). 

The binding energies in this figure are calculated including only the Coulombic interaction. 

The resonant formation rate of the complex in reaction (2) is very sensitive to the energy of the 
(1, 1) state of dtp . Small corrections due to relativistic, QED, hyperfine and other effects have 
been calculated. See, for example, refs. 29, 2 and 31. 

Reaction (2) is incomplete without an indication of what results from the formation of the 
[( dtp)ndee } complex. With inclusion of the most important decay products, it becomes 

tp + D 2 \{dtp)\\dee\ A user > [(dtp) Jv de] + 4- e, (3) 

back decay decay 


where ( J, v) = (0, 1), (2, 0), (1, 0) or (0, 0). If the total angular momentum of the tp + D 2 is 
taken to be zero for simplicity, the cross section a r (E) for this resonant process is determined 
by the Breit-Wigner formula (ref. 32). 


( p\ — — r e r a 

rK ' fc2 (£-£ r ) 2 + i ( r e + r a ) 2 ’ 


(4) 
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where 

E = the energy of the system, 

E r ~ the energy of the resonant state, 
k = wave number of the relative motion of tp and D 2 , 
r e = partial width for back decay, 

T 0 = partial width for Auger decay. 

p 

It follows from this that the deexcitation rate of [( dtpn)dee ] — It is very important that 
this is sufficiently fast that the [{dtp)\\dee\ does not simply decay back into tp + D 2 (ref. 33). 

Bhatia et al. (ref. 34) carried out a calculation of the deexcitation rate of the (1, 1) state of dtp 
in [(dtp)udee] into all states of lower energy. In [{dtp)\\dee\, the d , t and p are close together. 
To a good approximation, they can be looked upon as forming a nucleus of charge +1 in a 
molecular complex, similar to D 2 , but with one d nucleus replaced by a fictitious particle with 
mass equal to that of dtp. 

To be able to carry out their calculations, Bhatia et al. simplified this to a dtp in its (1, 1) state 
acting as the nucleus for a hydrogen- like atom in its ground state. They took the initial and 
final wave functions 1 fa and to be of the form 

il>i = / 0»(fi,r 2 )^i s (r), 
t \>f =^/(ri,r 2 )F fc (r), 

where r\ and r 2 are the position vectors of the d and p, respectively, with respect to the t as 
origin and 

v = position vector of the electron with respect to the center of mass of dtp as origin. 

^i(r i»r 2 ) and i/7(r l5 r 2 ) are the initial and final state wave functions for the dtp . They were 
calculated using Hylleraas-type basis functions in a similar way to the calculation described in 
ref. 19. 

jpi s (r) is the Is initial state of the electron and Ffc(r) is its final continuum state. 

The Auger decay of the dtp in its (1, 1) state is brought about by the finite size of the dtp . If 
the d, t and p were all at the same point, the system would indeed be equivalent to a hydrogen 
atom in its ground state. However, this is only approximately the case; a finite size correction 
must be included. This is of the form of a multipole expansion 

V = — 5 - • [(mi - l)r a + (m 2 + l)r 2 ] + . . . (5) 

r* 

where e is the charge on the proton, 

mi = 77 , 7712 — ”TT" (M = 771^ + 77l t -f 771 J 

M M 

and md, for example, is the mass of the deuteron. 

The V can be treated as a perturbation using Fermi’s Golden Rule. See, for example, ref. 35. 
The results Bhatia et al. obtained by including the first, dipole term in (5) are given in Table 
1. Note that the deexcitation rate to the (0, 1) state with zero angular momentum is much 
larger than the other deexcitation rates. 

The results obtained in some other calculations for the key (1, 1) — ► (0, 1) deexcitation rate are 
given in Table 2. The calculation in which I was involved (ref. 38) differs from the others in 
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that it allows for the molecular structure of the [(dt^)ndee\ complex. It can be seen that 
allowing for this does not significantly affect the deexcitation rate. 

The circumstances that led to this calculation are of interest. In December 1987, I visited the 
atomic physics group at Goddard Space Flight Center to give a seminar on calculations I was 
carrying out on positron-hydrogen- molecule scattering. Drachman talked to me about the 
calculations that he was carrying out with Bhatia and Chatterjee (ref. 34) on the deexcitation 
rate of the [{dti±)\\dee] complex. Knowing that I had come into positron physics from 
theoretical chemistry, he encouraged me to see if I could take the molecular nature of the 
complex into account when calculating the deexcitation rate. As I pointed out earlier, this was 
not allowed for in ref. 34. He gave me a translation that he had made of a report by two 
members of Ponomarev’s group, Men’shikov and Faifman, on the treatment of [(dtfi)ndee] as a 
molecular complex. 

I was sufficiently interested to put my graduate student, Lewis, onto the problem of how to 
take into account the molecular nature of [(dtfi)udee\ when calculating the deexcitation rate. 

It was not long before I came to appreciate the special features of muon catalyzed fusion that 
had attracted Drachman to it. 

Lewis and I collaborated with Hara at RIKEN in Japan to allow for the fact that [{dt^i)dee\ 
and [(dt(j>)de] + in reaction (3) should be treated as D 2 -like and Dj -like molecules, respectively, 
with one nucleus replaced by a particle of mass equal to that of dt[i (ref. 38). 

Quite soon I was in contact with Ponomarev. I visited him and his group at the Kurchatov 
Institute in Moscow in 1990 and gave a seminar on my work. The Cold War had ended and 
contacts with Russian scientists were multiplying to make up for the long period of separation. 

With my postdoc, Harston, I went on to collaborate with Faifman and Strizh on the 
calculation of the rate of the reaction (3) that leads to the formation of dtp in a state with zero 
angular momentum. Men’shikov and Faifman (ref. 32) considered that the inequalities 
T c <C T a e, where c is the kinetic energy of the relative motion of tfi and D 2 , held for T a and 
r € in equation (4). Under these conditions, cr r (E) is effectively the cross section for the 
formation of the complex (ref. 39). Also, 


a r n-jp-Tt5(E-Er). ( 6 ) 

T e is determined by the matrix element (x/|F|xi), where Xi and Xf 3X6 the initial state _ 
tfi + D 2 wave function and the wave function for the [(dtfi)udee\ complex, respectively. V is 
the potential that brings about the formation of the complex, i.e. the interaction between the 
f/i and D 2 . Men’shikov and Faifman showed that it could more conveniently be taken to be the 
analogue of V in equation (5), allowing for the molecular structure of the complex. 

By skillful use of mathematics (ref. 40), Faifman et al. (ref. 41) were able to evaluate (Xf\V\Xi)) 
with V the leading dipole term in its multipole expansion. During visits to Nottingham by 
Faifman and Strizh, this was extended to include the next term which is the quadrupole term 
(ref. 42). It was found that including the quadrupole term reduced the magnitude of the peak 
reaction rates by between 20% and 30%. 

In fact, this treatment does not allow fully for the molecular nature of the resonant process in 
which dtfjL is formed in a state with J = 0 in which fusion can take place rapidly. Cohen at Los 
Alamos, who has made many contributions to the theory of muon catalyzed fusion over the 
years, encouraged me to apply the methods of quantum reactive scattering to reaction (3). 

This reaction has some similarities with the chemical reaction 

H + D 2 — > H + D 2 or HD + D. (7) 
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This and other chemical reactions have been extensively studied by these methods. See, for 
example, ref. 43. 

However, comparison of reactions (3) and (7) shows that (3) has a special feature not present 
in (7), namely the Auger decay process which leads to the loss of the electron. Reaction (7) 
corresponds quite closely to the much slower side reaction 

t/i + D 2 ±=+ {(dtp) nee] — ► [(dtp)j v e] + D (8) 

where (J,v) = (0,1), (2,0), (1,0) or (0,0). 

This reaction proceeds through the same resonances as reaction (3). In this case, the resonant 
complex decays forward into the muonic molecule \(dtp)j v e\ and a D atom or backward into 
tfj, + D 2 . 

One of the methods that has been applied to reaction (7) is the method of Pack and Parker 
(ref. 44). They use adiabatically adjusting, principal axes hyperspherical coordinates (APH) in 
their calculations. These are elegant coordinates that transform smoothly between different 
channels such as H + D 2 and HD + D. Together with Pack, my postdoc Zeman and I applied 
this method to reaction (8) (refs. 45-47). 

Unfortunately, there was no easy way of including the Auger decay channel directly in our 
treatment. However, as pointed out by Men’shikov and Faifman (ref. 32), the coupling between 
the resonant channels is small as the lifetime of the resonant complex is much longer than the 
time the complex takes to complete a vibration. Thus the various decay processes operate 
essentially independently. This made it possible for us to obtain what we expect to be accurate 
values for T e and r n , the partial widths for back decay and for the decay of the complex into 
[(dtp) j v e\ + D, by analysing our results for the cross section for the resonant reaction (8) using 
the Breit-Wigner formula. 

Owing to the complexity of the calculation, we were only able to consider scattering states with 
zero total angular momentum. As the dtp is in its (1,1) state in the complex, this meant that 
the fictitious diatomic molecule in the complex, with one nucleus dtp, had to be in a state with 
angular momentum equal to 1. Care had to be taken to make sure that all important 
corrections to the Born-Oppenheimer potential were included in calculating the potential 
energy surface for the reaction. 

An interesting picture of the surface function eigenvalues of the hyperradius p is given in 
Figures 3 and 4 in ref. 46. In Figure 3, the bottom four curves correspond asymptotically to the 
open channel (0, 1), (2,0), (1,0) and (0,0) states of [( dtp)j v e ]. Figure 4 highlights the curve 
that corresponds to the closed channel (1,1) state asymptotically among curves that 
correspond in this region to the various D 2 states. 

It was found, as expected, that the reaction proceeded only through the resonances. The 
resonant states were found to be the vibrational states of the complex with vibrational 
quantum number v c = 3 and 4. The center of the v c = 2 resonance was calculated to be just 
slightly below threshold. 

Somewhat to our surprise, we found that our calculated value of the back decay rate (ref. 47) 
was much larger than the value obtained for it by Lane (ref. 33) and comparable with the 
calculated values of the Auger decay rate in Table 2. Further work is necessary to resolve this 
discrepancy. 

It was a great pleasure for me to be able to give Drachman a copy of this treatment of the 
resonant process that leads to fusion (ref. 46). It showed what had resulted from his 
encouragement to me to consider the molecular nature of the complex in reaction (3) . 
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MUON CATALYZED FUSION AS A SOURCE OF ENERGY 

I have not said anything so far about muon catalyzed fusion as a possible source of energy. To 
consider this, we must look at the full muon catalyzed fusion cycle. This begins with the 
projection of a muon into a mixture of deuterium and tritium at room temperature. The 
sequence of events is then as shown in Figure 2. 

What I have not considered up until now is how the cycle continues after fusion takes place. 
The muon is a lepton and is thus unaffected by the strong interaction that brings about the d-t 
fusion 

d 4 - 1 — ► a + n. (9) 

However, the muon is negatively charged and may combine with the a particle to form a bound 
state. This process is referred to as ‘sticking’. If the muon remains bound in this way, it is lost 
to the cycle. 

However, the afi is emitted with a kinetic energy of 3.5 MeV (ref. 30). As it slows down in the 
medium, it undergoes many collisions which may result in the stripping of the muon from the 
a particle. In addition, the muon can become available through transfer reactions that result in 
tfi or and a free a particle. 

We can define an effective sticking probability, w s , by 

w a - ti;°(l - R), (10) 


where w° is the initial sticking probability and R is the muon reactivation coefficient that 
measures the extent to which stripping occurs, has been calculated by several methods: 
variational (refs. 49-51), adiabatic (ref. 52), adiabatic hyperspherical (ref. 53) and Monte Carlo 
(ref. 54). Detailed calculations have been carried out of R by Cohen (ref. 55), Markushin 
(ref. 56) and Rafelski et al. (ref. 57). 

If we consider the rate values given in Figure 2, the muonic atom formation rate A a is much 
larger than either A dt , the rate of transfer of a muon from d/i in its ground Is state to tfx or 
Xdti . n the resonant formation rate of dt\x. Consequently, the rate of the muon catalysis cycle A c 
is practically independent of A a . For the same reason, it is practically independent of the dtfjL 
deexcitation rate (see Table 1) and the fusion rate. 

A c is then given by 

(u) 

where Cd and c* are the concentrations of deuterium and tritium, respectively, in the mixture 
(cd +Q = 1). See refs. 29, 48 and 58. qu is the probability that the dfi reaches its ground state. 
The rate of transfer of the muon from d to t is much slower for d/i in this state. The first term 
on the right-hand side of equation (11) is the time the muon spends as dfi waiting to transfer 
to t and the second term is the time it spends as tfi waiting to form the dtfi molecule. 

If sticking did not occur the average number X c of cycles catalyzed by one muon would be 

A c / Aq, where — = 2.2 x 10 -6 sec is the mean lifetime of the muon. Sticking modifies it to 
Ao 


*c = 


1 



( 12 ) 


Theory and experiment indicate that w s « 0.5% (refs. 48, 58, 59). If w s — 0.43%, as in Figure 2, 
it follows from equation (12) that X c cannot be larger than 232. 
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With the estimated cost of ~ 8 GeV to produce a muon (ref. 60) and a maximum energy 
release of 17.6 MeV per cycle, break even would be achieved if X c « 450. A considerably 
higher X c value, at least 900, would be required for muon catalyzed fusion to be a 
commercially viable source of energy. X c = 150 is currently attainable. Possible ways of 
making X c as high as 900 by increasing the density of the deuterium and tritium mixture and 
decreasing by increasing R are outlined in ref. 61. See equation (10). 

14.1 MeV of the 17.6 MeV generated by the d-t fusion reaction (9) is carried away by the 
neutron. The present maximum X c value of ~ 150 would be adequate for the production of 
energy by using neutrons obtained in this way to convert 238 U into plutonium fuel for 
conventional fission reactors (refs. 29, 60, 30,62). Unfortunately, this method of energy 
production would not be as ‘clean’ as using the muon catalyzed fusion cycle directly. 

The greatly expanded interest in muon catalyzed fusion in the 1980s was followed by a gradual 
decline in interest with the realisation that it was unlikely to be a viable source of energy. At 
the beginning of the 90s, the Advanced Energy Projects (AEP) Division of the US Department 
of Energy withdrew funding from experiments on muon catalyzed fusion at Los Alamos. The 
AEP Division considered that the research had been successful as it had shown that muon 
catalyzed fusion was not a viable source of energy. 1 

Experimental work continued at the Paul Scherrer Institute (PSI) in Switzerland (ref. 63, 64) 
and at Dubna and elsewhere in Russia (ref. 65). In the late 90s and early 2000s, a considerable 
impetus was given to experimental work on muon catalyzed fusion by Nagamine and his 
Japanese group from RIKEN who carried out experiments using the pulsed muon beam at the 
Rutherford Appleton Laboratory (RAL) in the UK (refs. 59, 61, 66) and also by the 
international collaboration based on Marshall’s group at TRIUMF in Canada (refs. 67, 68). 
This essentially brings us to the present day. 

The reader who wishes to obtain an overall view of muon catalyzed fusion should choose from 
the review articles in refs. 2, 10,29,30,48,58,63,69 and 70. The most recent, detailed overview 
is given in the Proceedings of the International RIKEN Conference on Muon Catalyzed Fusion 
and Related Atoms (MuCFOl) at Shimoda in Japan in April, 2001 (ref. 71). 

CONCLUSION 

Today with the maximum number of fusions attained per cycle still only about 150, it is 
difficult to argue that muon catalyzed fusion should be considered seriously as a possible 
energy source. However, if Nagamine (ref. 61) or anyone else were able to increase this number 
to ~ 900, the situation would change, but this seems unlikely. 

This is unfortunate but over the years the study of muon catalyzed fusion has given many 
physicists, not least Dick Drachman and myself, many moments of real pleasure. Nobody 
knows why physical phenomena can be described using mathematics (ref. 72) but to apply it 
successfully to an elegant physical problem is an experience to be savoured. 

I am very grateful to Dick Drachman for introducing me to muon catalyzed fusion and for 
many valuable discussions with him over the years about it, positron physics and other 
subjects. It has always been a great pleasure for me to visit Goddard Space Flight Center and 
discuss my work with him, Aaron Temkin and Anand Bhatia. I wish Dick and Aaron all the 
best in their retirement. 


1 James Cohen, private communication, 2005. 
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Table 1 Calculated values for the deexcitation rates (ref. 34) 


Transition 

— 
Energy difference 

A E (eV) 

Deexcitation 
rate (sec“l) 

11 — ►01 

34.1745 

6.44 x 10 n 

o 

o 

T 

O 

j-H 

86.6682 

3.31 x 10 10 

11 20 

101.88 

8.56 x 10 10 

o 

rH 

T 

129.93 

5.12 x 10 10 

01 -4 10 

197.6370 

3.32 x 10 10 

11 ^00 

318.4797 

1.57 x 10 9 
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Table 2 Other calculations of the (1,1)-+ (0, 1) deexcitation rate 


Calculation 

Deexcitation rate 
11 — ► 01 (sec -1 ) 

Vinitsky, Ponomarev 
and Faifman (ref. 36) 

11.4 x 10 11 

Scrinzi and 
Szalewicz (ref. 37) 

10.20 x 10 u 

Armour, Lewis 
and Hara (ref. 38) 

8.63 x 10 u 

Bhatia, Drachman 
and Chatterjee (ref. 34) 

6.44 x 10 11 
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Figure 1 Energy level diagram and dipole Auger (or internal conversion) transitions f< 
dt/i molecular ion (ref. 30). Binding energies for J = 0 and J = 1 are from ref. 25 and for 
ref 26. Units are eV. 
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Figure 2 The principal muon catalysis fusion cycle in a deuterium and tritium mixture 
(ref. 29). Side chains involving d-d and t-t fusion are not shown. A a = muonic atom formation 
rate ^4x 10 12 a" 1 ; A dt — muon transfer rate from (dfi) i 5 to (tfj)u w 3 x 10 8 s' 1 ; A dtp ~ 
resonant formation rate of dtfi w4x 10 8 sec' 1 ; = effective sticking probability « 0.43%; 

A c = cycle rate. Values (ref. 48) are for T — 300K and liquid hydrogen density (4.25 x 10 28 
atoms m~ 3 ). 
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